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Abstract 

This  is  a  final  report  for  EOARD  grant  number  FA8655-05- 1-3054,  “Extension  of 
electromagnetic  code  for  phased  array  analysis  and  design”.  The  name  of  the  electromagnetic 
code  is  PB-FDTD.  The  unique  property  of  the  code  is  that  it  can  perform  unit  cell  analysis  of 
periodic  structures  (such  as  array  antennas)  in  the  time  domain.  This  means  that  a  broadband 
response  is  obtained  in  a  single  computation.  The  antenna  element  feed  is  a  key  part  of  phased 
array  antennas.  An  accurate  model  for  the  feed  is  therefore  important.  Some  type  of  feeds  can 
already  be  modeled  in  the  PB-FDTD  code.  Waveguide  ports  are  very  common  in  array  antennas. 
Such  feeds  have  however  not  been  implemented  in  the  PB-FDTD  code.  The  main  objective  of  the 
work  was  therefore  to  implement  waveguide  port  in  the  PB-FDTD  code.  This  has  been 
accomplished.  Another  objective  was  to  simplify  the  analysis  of  arrays  with  circular  polarized 
antenna  element.  This  and  more  have  been  done.  The  code  can  now  calculate  both  the 
polarization  and  the  transmission  coefficient  to  all  propagating  Floquet  modes.  The  GUI  for  the 
code  has  also  been  modified  in  order  to  handle  the  new  functions. 


Background 

Many  electromagnetic  problems  have  a  periodic  structure.  Examples  are  phased  array  antennas, 
frequency  selective  surfaces  (important  for  stealth  application),  meta-materials  (currently  a  hot 
research  area)  and  electromagnetic  bandgap  materials.  Electromagnetic  analysis  and  design  of 
such  structures  are  greatly  simplified  by  the  use  of  a  mathematical  theorem  called  Floquets 
theorem  [1].  With  Floquets  theorem  only  one  unit  cell  in  the  periodic  structure  needs  to  be 
modeled.  The  electromagnetic  interaction  to  all  other  unit  cells  is  automatically  taken  into 
account.  Previously,  unit  cell  analysis  based  on  Floquets  theorem  was  always  performed  in  the 
frequency  domain.  This  means  that  only  one  frequency  point  at  the  time  can  be  analyzed.  In 
many  cases,  a  broadband  response  is  required  (for  example  radar  cross  section  calculations  and 
broadband  array  design).  A  time  domain  version  of  Floquets  theorem  is  therefore  highly  desirable 
since  that  would  generate  data  for  all  frequencies  in  a  single  computation.  A  time  domain  version 
of  Floquets  theorem  was  developed  by  Henrik  Holter  during  his  PhD-studies  1997-2000  [2].  The 
new  method  has  been  implemented  in  the  very  popular  Finite-Difference  Time-Domain  method 
(FDTD)  for  solving  Maxwells  equations.  The  name  of  the  computer  code  is  PB-FDTD  (Periodic 
Boundary  FDTD).  The  use  of  FDTD  means  that  complicated  structures  and  materials  can  be 
modeled.  Some  years  ago,  EOARD  supported  the  development  of  a  graphical  user  interface  etc 
(GUI)  for  the  PB-FDTD  code  under  contract  no.  F61 775-00- WE037.  The  PB-FDTD  code  has 
after  that  continuously  been  improved.  It  can  now  handle  all  the  periodic  structures  mentioned 
above.  The  code  is  numerically  highly  efficient  and  has  proved  quite  useful  at  the  AFRL, 
Hanscom,  MA  and  other  places. 
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Objective 

The  PB-FDTD  code  has  primarily  been  used  for  broadband  phased  array  antenna  design  (several 
scientific  papers  have  been  published  in  periodicals  where  the  PB-FDTD  code  has  been  used). 

The  feed  is  a  key  part  of  array  antennas  and  therefore  needs  to  be  included  in  the  computer 
model.  Some  type  of  feeds  (coaxial  cables,  micro-strip  lines  and  striplines)  can  already  be 
modeled  in  the  PB-FDTD  code.  However  one  of  the  most  common  feeds  in  array  antennas, 
waveguide  ports,  cannot  be  modeled.  The  main  proposed  research  effort  is  therefore  to 
implement  waveguide  ports  in  the  computer  code.  That  would  make  the  code  still  more  useful. 
Waveguides  are  complicated  since  the  electromagnetic  field  inside  the  guide  is  described  by  a 
sum  of  forward  and  backward  traveling  propagating  and  non-propagating  modes.  The  shape  and 
frequency  dependent  propagation  constant  of  the  modes  depend  on  the  waveguide  cross-section 
and  the  material  inside  the  guide  and  must  in  general  be  calculated  numerically.  The  waveguide 
must  be  terminated  in  a  load  that  absorbs  all  backward  traveling  propagating  and  non-propagating 
modes.  Further,  the  use  of  waveguide  modes  together  with  unit  cell  technique  caused  an  intricate 
problem  that  had  to  be  solved  (described  later  in  the  text). 

An  additional  part  of  the  proposal  is  to  simplify  the  analysis  of  arrays  with  circular  polarized 
antenna  elements.  For  characterization  of  such  elements  with  the  current  code,  several 
computations  with  subsequent  post-processing  are  necessary.  The  proposed  research  effort  is  to 
simplify  this. 

Statement  of  work 

The  proposed  services  in  the  application  to  EOARD  are  described  in  this  section. 

The  waveguide  implementation  consists  of  the  following  parts: 

•  The  cross-section  of  the  waveguide  can  be  circular  or  rectangular.  These  are  the  most 
common  shapes.  The  code  will  almost  certainly  handle  almost  arbitrary  cross-sections,  but 
that  can  not  be  guaranteed.  (It  turned  out  that  the  code  can  handle  arbitrary  cross-sections). 

•  The  waveguide  can  be  filled  with  a  lossless  dielectric  material. 

•  The  waveguide  feed  is  directed  in  the  array  normal  direction. 

•  The  number  of  waveguide  modes  can  be  specified  by  the  user.  The  code  automatically 
calculates  the  modes  and  sorts  them  after  the  value  of  the  associated  transversal  wavenumber. 

•  The  user  specifies  which  of  the  modes  that  excite  the  guide. 

•  The  coupling  between  the  exited  mode  and  the  other  waveguide  modes  is  automatically 
calculated. 

•  For  degenerated  modes,  the  user  can  specify  any  linear  combinations  of  such  modes  in  order 
to  control  the  polarization. 

•  A  unit  cell  in  the  array  contains  one  waveguide  feed.  More  than  one  waveguide  feed  will 
probably  bee  implemented,  but  that  can  not  be  guaranteed. 

•  Examples  with  waveguide  feeds  is  included  in  the  manual  for  the  code 

•  The  array  lattice  is  rectangular.  Attempts  will  be  done  to  implement  waveguide  feeds  for 
triangular  lattices  also. 
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•  The  waveguides  can  also  be  used  for  scattering  (RCS)  calculations.  A  plane  wave  incident  on 
the  array  is  used  for  that. 

The  circular  polarized  element  analysis  part  consists  of  the  following  parts: 

•  Simplification  of  the  calculation  of  the  scattering  parameters  between  the  antenna  element 
and  the  zero  order  TE-  and  TM-mode  Floquet  modes.  It  is  accomplished  in  the  following 
way:  the  antenna  is  excited  and  the  coupling  to  the  zero  order  TE-  and  TM-modes  is 
calculated.  (It  turned  out  that  the  coupling  to  all,  not  only  the  zeros  order,  propagating  Floquet 
modes  is  calculated). 

•  An  example  in  the  manual. 

In  addition  to  the  above,  the  graphical  user  interface  (GUI)  will  be  updated  so  it  can  handle  the 
new  functionalities. 

Deliverables 

•  Computer  code,  inclusive  source  code  (FORTRAN  and  MATLAB). 

•  Updated  manual.  New  examples. 

Implementation  of  waveguide  ports 

In  this  section,  the  implementation  of  waveguide  ports  in  the  PB-FDTD  code  is  described  in 
detail. 

There  are  several  ways  to  implement  waveguide  ports  in  the  FDTD  method.  Some  research  has 
therefore  been  carried  out  in  order  to  find  the  most  accurate  and  usable  method.  The  workload  for 
implementing  the  method  has  not  been  a  constraint.  The  chosen  method  is  described  in  [3].  It  is  a 
mod-matching  technique. 

Since  a  mode  matching  technique  is  used,  evanescent  modes  can  be  included  in  the  expansion  of 
the  field  in  the  waveguide.  This  means  that  the  port  does  not  need  to  be  very  long  in  order  to 
correctly  handle  evanescent  modes. 

There  are  two  steps  in  order  to  implement  the  mod-matching  technique.  First  the  waveguide 
modes,  transversal  wave-numbers  and  the  time  domain  impulse  response  of  each  mode  must  be 
determined.  Second,  the  waveguide  modes  together  with  the  impulse  responses  are  used  to 
truncate  and  excite  the  waveguide  in  the  FDTD  method.  These  two  steps  are  described  below. 

Waveguide  modes 

Since  the  waveguide  cross  section  can  have  an  almost  arbitrary  shape,  the  waveguide  modes  must 
be  determined  numerically  (as  a  matter  of  fact,  even  if  the  waveguide  where  rectangular,  using 
the  known  theoretical  modes  and  wave-numbers  would  result  in  inaccurate  results.  This  has  to  do 
with  the  mesh  in  FDTD).  In  order  to  find  the  modes,  Helmholtz’s  two-dimensional  equation  for 
the  longitudinal  part  of  the  electric  or  magnetic  field  or  Laplace  equation  for  the  potential  is 
solved  together  with  suitable  boundary  conditions. 


Page  3 


25  October,  2005 


(v2  +  k2  )hz  =  0,  TE  -  modes 

(v2  +  k2r  )EZ  =  0,  TM  -  modes  (1) 

V2F  =  0,  TEM-modes 
The  boundary  conditions  are 

=  0,  TE  -  modes 
dn 

Ez  =  0,  TM  -  modes  (2) 

F  constant  on  each  conductor,  different  on  different  conductors 

Solving  (1)  on  the  FDTD-mesh  results  in  a  matrix.  For  TE-  and  TM-modes,  the  eigenvalues  and 
eigenvectors  of  this  matrix  correspond  to  the  transversal  wave-numbers  squared  and  the 
longitudinal  magnetic  or  electric  field.  For  TEM-modes,  the  potential  in  the  waveguide  cross 
section  is  obtained. 

The  transversal  electric  field  is  needed  for  the  mode  matching  technique.  For  TE-  and  TM-modes 
the  transversal  field  is  obtained  via  Maxwell’s  equations.  For  TEM-modes,  the  electric  field  is 
obtained  from  the  gradient  of  the  potential 

£r=-VF  (3) 

Equation  (1)  is  solved  on  the  same  mesh  as  the  FDTD-mesh.  That  guarantees  that  the  accuracy 
will  be  high.  Observe  that  (1)  is  a  frequency  domain  equation  for  the  modes.  The  modes  are 
however  the  same  in  the  time  domain. 

Construction  of  the  matrix  above  for  arbitrary  waveguide  cross-sections  and  taking  into  account 
the  boundary  condition  turned  out  to  be  a  complex  task.  A  MATLAB  and  a  FORTRAN  program 
have  been  developed  for  this  part. 

The  time  domain  impulse  response  of  the  waveguide  modes  is  also  required  for  the  mode 
matching  technique.  In  [3],  it  is  showed  how  the  impulse  responses  can  be  calculated  numerically 
by  a  recursive  formula.  This  method  has  been  used  here. 

Implementation  in  FDTD 

In  this  section,  the  implementation  of  the  mode  matching  technique  in  the  FDTD  method  is 
described  step  by  step. 


FDTD  mesh 


Port  position 

Figure  1.  Waveguide  example.  Side  view. 
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The  electric  field  in  the  waveguide  cross-section  at  z=0  in  Figure  1  is  going  to  be  computed  from 
the  known  electric  field  in  the  cross-section  at  z=Az  and  the  incident  signal  at  z=0.  Applying  the 
computed  field  at  z=0  will  both  truncate  and  excite  the  waveguide.  The  truncation  is  in  theory 
reflection  free.  Examples  later  in  this  report  show  that  the  reflection  is  very  small. 

Assume  that  the  FDTD  cell  size  in  the  z-direction  is  given  by  A z  and  the  discrete  time  4  in  the 
FDTD  method  is  given  by 

fr  =  kAt ,  (4) 


where  At  is  the  FDTD  time  step. 

The  total  transversal  electric  field  at  z  =  Az  can  be  expressed  as  an  expansion  of  N 
orthonormalized  waveguide  modes  (evanescent  modes  can  be  included) 

Et  (A Z,  kAt )  =Yjan  (kA&n  (5) 

n=l 

Er(Az,kAt)  is  a  known  quantity  (from  the  FDTD  method).  The  voltage  an  is  determined  from 
orthogonality 

an (kAt)  =  ^Et (Az, kAt ) •  ends,  n-\..N  (6) 

s 


To  obtain  the  backward  voltages,  a.„,  the  incident  voltage  needs  to  be  subtracted  from  the  a„. 

Assume  that  the  incident  signal  at  z=0  is  given  by  s(t)  and  that  mode  m  is  the  source  (incident) 
waveguide  mode. 


The  incident  signal  at  z=Az  is  given  by  the  convolution  between  s(t)  and  the  waveguide  impulse 
response  for  mode  m. 


k- 1 


$&z  (kAt)  =  Pm(k)*  s(kAt)  =  £  s(qAt)Pm  (k-q) 

9=1 

The  impulse  response  Pm  is  calculated  according  to  [3]. 
The  backward  voltages  at  z=Az  is  given  by 

[an(kAt)-s^(kAt\  n  =  m 


a-SkAt)  = 


an(kAt),n  *  m 


,  n  =  \..N 


(7) 


(8) 


The  backward  voltages  at  z=0  is  given  by  the  convolution  between  the  a.„  and  the  impulse 
response. 

<„° (kAt)  =  Pn (k) * a_n (kAt)  =  £ a_n (qAt)Pn (k-q),  n  =  \..N  (9) 

9=1 

Next,  the  voltages  of  the  total  field  at  z=0  is  calculated  by  adding  the  incident  signal  s(t)  to  (9) 
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z= 0 


(Mr)  = 


\  az_  „°  (Mr) + s(kAt),  n-m 
a!"„°  (Mr)  n^m 


n  =  l..N 


(10) 


The  total  transversal  field  at  z=0  is  finally  given  as  a  superposition  of  waveguide  modes 

£r(z=o,M>)=y<7f"(ta/)>,  (ii) 

«= i 

This  field  is  applied  at  the  plane  z=0  in  the  FDTD-code.  As  a  result,  backward  traveling  modes 
will  be  absorbed  at  z=0  (if  they  are  included  in  the  N  modes).  At  the  same  time,  the  excitation 
signal  s(t )  is  applied  at  z=0. 

The  above  steps  are  repeated  for  each  time  step  in  the  FDTD  method.  The  method  requires  that 
all  a_  n  (Mr)  and  a*=0(kAt)  is  saved  for  all  FDTD  time  steps  k.  Memory  is  also  required  to  save 

the  impulse  response  for  each  waveguide  mode  and  the  excitation  signal.  The  memory 
requirement  for  all  this  is  however  small  compared  with  the  memory  requirement  for  the  FDTD 
mesh. 


Some  remarks  on  the  shape  of  the  excitation  signal  s(t ):  The  spectrum  of  the  pulse  should  have  an 
energy  distribution  that  is  adapted  to  the  desired  frequency  range.  It  is  very  common  to  use  a 
Gaussian  pulse  in  the  FDTD  method.  A  waveguide  however  has  a  cut-off  frequency.  The  pulse 
should  therefore  contain  little  energy  below  and  near  the  cut-off  frequency.  Otherwise,  very  long 
simulation  times  caused  by  slowly  propagating  signals  just  above  the  cut-off  frequency  may 
occur.  In  the  literature,  people  have  suggested  to  use  a  sine-pulse  in  the  time  domain  (square 
pulse  in  the  frequency  domain)  for  waveguide  excitation  in  the  FDTD  method.  This  has  been 
tried  here.  It  was  however  found  that  such  pulses  caused  some  problems.  One  problem  is  that  the 
shape  of  the  sine -pulse  in  the  time  domain  decay  slowly  (1  /t).  Since  a  truncated  sine-pulse  must 
be  used  in  the  FDTD  method,  the  Fourier  transform  will  not  be  a  perfect  square.  And  it  will 
therefore  be  energy  below  the  cut-off  frequency.  Second,  using  a  sine-pulse  with  the  unit-cell 
technique  in  PB-FDTD  sometimes  resulted  in  numerical  instabilities.  The  reason  for  that  is  again 
that  a  truncated  sine-pulse  must  be  used  (the  sine-pulse  is  not  very  small  when  the  computation 
starts).  A  Gaussian  pulse  has  an  exponential  decay  in  both  the  time-  and  the  frequency  domain. 
This  is  much  better  for  the  unit  cell  technique.  A  Gaussian  pulse  is  therefore  used  in  PB-FDTD 
for  excitation  of  waveguides.  When  the  user  specifies  the  frequency  range  in  PB-FDTD,  the  start 
frequency  should  not  be  below  the  cut-off  frequency. 


Mode  matching  and  the  PB-FDTD  unit  cell  technique 

The  mode  matching  technique  require  that  the  field  in  the  cross-section  of  the  waveguide  at  z=Az 
is  available.  This  result  in  a  difficulty  for  the  unit  cell  technique  in  PB-FDTD.  The  unit  cell 
technique,  described  in  [4],  is  based  on  a  unit  cell  that  during  the  computation  is  moved  over  the 
periodic  structure  (for  example  an  array  antenna)  with  at  least  the  speed  of  light.  This  is 
illustrated  in  Figure  2  where  each  unit  cell  has  a  circular  waveguide.  The  dark  unit  cell  in  Figure 
1  is  the  unit  cell  modeled  in  FDTD.  It  is  that  unit  cell  that  is  moved  over  the  surface  during  the 
computation. 
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Unit  cell  p-1  p  p+1  p+2 

number 


Figure  2.  Unit  cell  technique.  A  moving  unit  cell  is  used.  The  circle  in  each  unit  cell  is  a  waveguide  port. 

From  equation  (6)  in  the  mode  matching  technique,  it  is  obvious  that  the  transversal  field  in  the 
waveguide  at  z~Az  is  needed.  The  problem  is  that  according  to  Figure  2,  only  a  part  of  the 
waveguide  may  be  inside  the  moving  unit  cell  (it  is  not  possible  just  to  combine  the  two  partial 
visible  waveguides  in  the  moving  unit  cell  in  Figure  2  since  there  is  a  time  shift  between  them). 
The  total  transversal  field  in  the  waveguide  cross  section  at  z=Az  is  therefore  not  available  at  all 
FDTD  time  steps. 

A  non-approximate  solution  for  this  intricate  dilemma  was  however  found.  The  solution  to  the 
problem  is  first  to  use  two  instead  of  one  moving  unit  cell,  see  Figure  3.  That  guarantees  that  a 
whole  waveguide  is  visible  for  a  sufficient  long  time.  Second,  the  partially  visible  waveguide  in 
unit  cell  p+1  in  Figure  3  still  constitute  a  problem.  The  total  transversal  field  at  z=Az  is  not 
available.  So  it  is  not  possible  to  calculate  the  a_n(kAt )  and  az„=0(kAt)  from  (6)-(10).  That  is 
however  not  needed  since  these  values  have  already  been  calculated  for  the  waveguide  in  unit 
cell  p  at  time  kAt  -  St .  The  time  St  is  related  to  the  scan  angle,  0O ,  the  element  spacing,  d,  and 
the  speed  of  light  cq. 


(11) 


This  solves  the  problem.  So  as  long  as  the  rightmost  waveguide  is  just  partially  visible  in  the 
moving  unit  cell,  a_n{kAt)  and  azn=0(kAt)  are  obtained  from  the  waveguide  to  the  left  of  this 
waveguide  at  an  earlier  time  step. 
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Figure  3.  Two  unit  cells  modeled  in  FDTD. 


Frequency  domain  data 

Frequency  domain  data  is  obtained  via  a  Fourier  transform  of  the  backward  traveling  voltages 
az_^(kAt)  and  the  in-signal  s(t). 

The  coupling  between  the  excited  waveguide  mode,  m,  and  all  other  waveguide  modes  is 
calculated  as 


VZ  v. 

Kv- 


,n  =  l..N, 


(12) 


where  Vn  are  the  frequency  domain  modal  voltages.  Vm  is  the  Fourier  transform  of  the  in-signal 
s(t).  Z„  and  Zm  are  the  wave-impedances  of  mode  n  and  m.  The  wave  impedance  is  calculated 
from  the  transversal  wave-number,  kT,  which  is  obtained  when  equation  (1)  is  solved. 
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TM 


COE 

kT 


Z 


TEM  ~~ 


CO 
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(13) 


Normalization  with  the  square  root  of  the  wave-impedances  in  (12)  results  in  generalized 
scattering  parameters. 


Examples 

Examples  with  waveguide  ports  have  been  added  to  the  manual  for  the  PB-FDTD  software.  One 
example  is  a  phased  array  with  ridged  waveguides  radiating  into  free  space.  Results  from  this 
example  are  shown  below.  The  fundamental  mode  is  shown  in  Figure  4  (all  modes  can  be  plotted 
in  the  PB-FDTD  code).  The  generalized  scattering  parameters  are  shown  in  Figure  5.  A 
comparison  between  PB-FDTD  and  Ansoft  HFSS  is  shown  in  Figure  6  and  Figure  7.  Ansofit 
HFSS  is  a  popular  commercial  frequency  domain  code  that  can  handle  periodic  structures.  Since 
it  is  a  frequency  domain  code,  it  must  redo  the  calculation  for  each  frequency  point  of  interest. 
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Port  1,  TE-mode  1.  Cut-off:  1.659748  GHz 


Figure  4.  Fundamental  mode  in  the  waveguide. 


Frequency  (GHz) 

Figure  5.  S-parameters.  Mode  1  is  the  source  mode.  Scan  angle:  45  degrees  in  the  H-plane.  The  vertical  lines 
indicate  when  higher  order  modes  in  the  waveguide  start  to  propagate. 
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Figure  6.  Magnitude  of  Sn.  Comparison  with  Ansoft  HFSS.  Scan  angle:  45  degrees  in  the  H-plane.  First 
grating  lobe  at  3.2  GHz. 


Figure  7.  Phase  of  Sn. 

The  next  example  is  an  air-filled  rectangular  waveguide  that  instead  of  radiating  into  free  space  is 
terminated  in  a  short.  The  inner  dimension  of  the  waveguide  is  48.75x22.5  mm.  The  length  of  the 
guide  is  25  mm.  The  frequency  range  has  been  set  to  1-10  GHz. 

6  TE-modes  and  2  TM-modes  are  used.  TE-mode  1  is  the  source  mode.  Below  are  the  calculated 
cut-off  frequencies  for  the  modes: 

Mode  1,  TE-mode:  1,  Cut-off:  3.156995  GHz 
Mode  2,  TE-mode:  2,  Cut-off:  6.308597  GHz 
Mode  3,  TE-mode:  3,  Cut-off:  6.658208  GHz 
Mode  4,  TE-mode:  4,  Cut-off:  7.368742  GHz 
Mode  5,  TE-mode:  5,  Cut-off:  9.172248  GHz 
Mode  6,  TE-mode:  6,  Cut-off:  9.449421  GHz 
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Mode  7,  TM-mode:  1,  Cut-off:  7.368742  GHz 
Mode  8,  TM-mode:  2,  Cut-off:  9.172248  GHz 

The  calculated  scattering  parameters  are  shown  in  Figure  8  -  Figure  10.  The  coupling  from  TE- 
mode  1  to  higher  order  modes  should  in  theory  be  zero.  Because  of  numerical  noise,  there  is 
always  a  finite  coupling.  As  can  be  seen  in  Figure  8,  the  coupling  is  very  small  (<  -250  dB).  The 
vertical  lines  in  the  figures  indicate  where  higher  order  modes  start  to  propagate  in  the  waveguide 
(the  cut-off  frequencies  above). 


Port  1 


Figure  8.  Sn  -  Sg2. 

Sn  is  shown  in  Figure  9  and  Figure  10.  Sn  should  in  theory  be  equal  to  0  dB  (total  reflection 
because  of  the  short.  As  can  be  seen  in  the  figures,  Sn  is  very  close  to  0  dB.  In  Figure  9,  the 
Return  Loss  does  not  go  directly  to  zero  at  the  cut-off  frequency  (3.16  GHz)  of  the  first  TE-mode. 
This  is  almost  certainly  because  of  very  slowly  propagating  signals  at  frequencies  close  to  the 
cut-off  frequency. 


Port  1 


Figure  9.  Sn. 
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Figure  10.  Sn.  Observe  the  scale  on  the  y-axis. 

Third  example:  calculated  return  loss  for  a  single  polarized  Vivaldi-array  scanned  to  30  degrees 
in  the  H-plane  is  shown  in  Figure  11.  The  element  is  fed  by  a  stripline.  This  TEM-type  of 
transmission  line  could  already  be  modeled  with  another  (less  advanced,  more  lengthy)  technique 
in  PB-FDTD.  The  agreement  between  the  two  methods  in  Figure  1 1  is  very  good. 


Figure  11.  Return  loss  of  an  array  with  Vivaldi-elements.  Two  different  feed  techniques. 

The  Vivaldi-element  and  the  stripline  port  TEM-mode  is  shown  in  Figure  12  and  Figure  13. 
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Figure  12  A  tapered  slot  element  with  plated  through  vias  (one  of  the  metal  fins  has  been  removed  to  increase 
the  visibility).  The  element  is  fed  by  a  stripline. 
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Figure  13.  Calculated  stripline  port  TEM-mode.  Electric  field  in  the  cross-section  is  shown. 

Polarization  and  transmission  coefficient 

Two  parameter  of  considerable  interest  in  phased  array  antennas  is  the  polarization  and  the 
transmission  coefficient  between  the  antenna  feed  and  the  Floquet  modes.  Calculation  of  these 
parameters  has  been  implemented  in  the  PB-FDTD  code. 

The  original  objective,  as  mentioned  earlier  in  the  report,  was  to  simplify  the  analysis  of  arrays 
with  circular  polarized  antenna  element.  The  implementation  here  is  more  general  than  that. 

Polarization 

The  polarization  is  presented  in  the  form  of  magnitude  and  phase  of  the  electric  field  components 
for  all  propagating  Floquet  modes.  An  example  (dual  polarized  Vivaldi  phased  array  scanned  into 
the  diagonal  plane)  from  the  PB-FDTD  manual  is  shown  in  Figure  14  and  Figure  15  (no  higher 
order  propagating  modes). 
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Floquet  Mode  (0,  0) 


Frequency  (GHz) 


Figure  14.  Polarization.  Relative  magnitude  between  the  electric  field  components. 


Frequency  (GHz) 


Figure  15.  Polarization.  Phase  of  the  electric  field  components. 

Transmission  coefficient  calculation 

The  transmission  coefficient  between  an  antenna  feed  and  all  propagating  Floquet  modes  is 
calculated  according  to  the  following  equation 


Ji 

i,n  _  y 

“7? 


rm,n 


'Floquet  Y  Floquet 


Lx(y) 


Antenna  Antenna 


(14) 
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Zp£  is  the  wave  impedance  of  Floquet  mode  (m,n) 

Z Antenna  is  the  source  impedance.  For  point  sources:  equal  to  the  internal  resistance.  For 
waveguide  ports:  equal  to  the  wave-impedance  of  the  excited  mode. 

Vpi’"  is  the  amplitude  of  the  normalized  Floquet  mode  (m,n) 


V Antenna  is  the  in-signal  to  point  sources  or  waveguide  ports. 


Normalization  with  the  square  root  of  the  wave-impedances  in  (14)  results  in  generalized 
scattering  parameters. 

In  PB-FDTD,  structures  are  periodic  in  the  xy-plane.  The  z-direction  is  normal  to  the  surface  of 
the  periodic  structure.  The  transmission  coefficient  in  (14)  is  calculated  for  the  coupling  to  the  x- 
and  y-components  of  the  electric  field  of  the  Floquet  modes. 

For  loss-free  array  antennas,  experiment  calculations  have  shown  that  energy  conservation  is 
satisfied  (a  necessary  but  not  sufficient  condition  for  proving  that  everything  is  correct). 


lni!+slr"f+Zl7-r|!=l> 


(15) 


where  VA  is  the  array  active  reflection  coefficient.  The  summation  is  over  all  propagating  Floquet 
modes. 

An  example  (dual  polarized  Vivaldi  phased  array  scanned  into  the  diagonal  plane)  from  the  PB- 
FDTD  manual  is  shown  in  Figure  16  (no  higher  order  propagating  modes). 
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Figure  16.  Transmission  coefficient,  magnitude  and  phase. 


Conclusions 

More  than  all  of  the  proposed  services  in  the  statement  of  work  have  been  fulfilled.  It  was 
promised  that  rectangular  and  circular  waveguides  should  be  implemented.  The  code  can  actually 
handle  waveguides  with  arbitrary  cross-section.  It  can  also  handle  waveguides  with  TEM-modes 
(ex.  coaxial  cables  and  striplines). 

It  was  also  stated  that  analysis  of  circular  polarized  antenna  elements  should  be  simplified.  The 
implementation  in  the  code  is  more  general  than  that  (the  polarization  and  transmission 
coefficient  calculation). 

There  is  a  graphical  user  interface  (GUI)  for  the  code.  The  GUI  has  been  modified  in  order  to 
handle  the  new  functions. 
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